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ABSTRACT 

We present a new optical polarimetric catalog for the Small Magellanic Cloud (SMC). It contains 
a total of 7207 stars, located in the northeast (NE) and Wing sections of the SMC and part of the 
Magellanic Bridge. This new catalog is a significant improvement compared to previous polarimetric 
catalogs for the SMC. We used it to study the sky-projected interstellar magnetic field structure 
of the SMC. Three trends were observed for the ordered magnetic field direction at position angles 
(PAs) of (65° ± 10°), (115° ± 10°), and (150° ± 10°). Our results suggest the existence of an ordered 
magnetic field aligned with the Magellanic Bridge direction and SMC’s Bar in the NE region, which 
have PAs roughly at 115°.4 and 45°, respectively. However, the overall magnetic field structure is 
fairly complex. The trends at 115° and 150° may be correlated with the SMC’s bimodal structure, 
observed in Cepheids’ distances and HI velocities. We derived a value of B s y y = (0.947 ± 0.079) fiG 
for the ordered sky-projected magnetic field, and SB = (1.465±0.069) pG for the turbulent magnetic 
field. This estimate of H s k y is significantly larger (by a factor of ~ 10) than the line of sight field 
derived from Faraday rotation observations, suggesting that most of the ordered field component is on 
the plane of the sky. A turbulent magnetic field stronger than the ordered field agrees with observed 
estimates for other irregular and spiral galaxies. For the SMC the B s k y /5B ratio is closer to what is 
observed for our Galaxy than other irregular dwarf galaxies. 

Subject headings: galaxies: ISM, Magellanic Clouds, magnetic fields — techniques: polarimetric 


1. INTRODUCTION 

The Small Magellanic Cloud (SMC) and the Large 
Magellanic Cloud (LMC) are gas rich irregular galaxies 
and satellites of the Milky Way (MW). Given the prox¬ 
imity to the Galaxy and their high Galactic latitude, the 
Magellanic Clouds (MCs) are perfect laboratories for ex- 
tragalactic studies, since the light emitted from them is 
not much attenuated by the dust present in the Galactic 
disk. These objects are highly studied, motivated by the 
fact that, through understanding them, we can also bet¬ 
ter comprehend, for instance, galaxy evolution from high 
redshifts to now; the structure of irregular dwarf galax¬ 
ies; the interstellar medium (ISM) content at the external 
parts of huge spiral galaxies; star formation at low met- 
alicity environments; dwarf galaxies in interaction; and 
the relation of satellite galaxies with their hosts. 

The MCs are low metalicity galaxies 0. 5 and 


0.2ZO) for the LMC and SMC, respectively (Kurt & Du- 


four 1998). This peculiarity makes their ISM particularly 
different 


differences, 


i'om the one of our Galaxy. Among the major 


1 Stanimirovic et al.||2000 Bot et al. 2004 

Meixner et al. 

2010 

and the submm excess emission (1 

Bernard et al. 

2008 

Israel et al. 

2010: Bot et al. 2010; Planck Coilabo 2 

ration et ai.||2011 

). 


The Bar component is the main body of the SMC and 
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corresponds to its northern part, which has a position 
angle (PA) of ~ 45°. The Wind is located ~ 2° East of 


the Bar an d corresponds to its southern part ((van den 
B ergh|[2 007). 


The dynamical evolution of the SMC-LMC-MW sys¬ 
tem seems to have generated two interesting structures - 
the Magellanic Stream and the Magellanic Bridge. These 
structures are well illustrated in the HI co lumn density 
map (Figure 2) presented by Brims et al. (2005). The 
Magellanic Stream is a gas tail covering over 100°, also 
going through the Galactic south pole. One scenario says 
that this gas might have been pulled out of the SMC 
and LMC due to tidal interact ion between the MCs and 
the MW (Putman et al. 20031. On the other side, the 
Magellanic Bridge seems to have been generated due to 
the same eff ect, but between the SMC and the LMC 
(Bekki 2009). More recently, high precision proper mo 
tions from the Hubble Space Telescope showed that the 
MCs are moving 80 km/s f aster than earlier estimates 
(Kallivayalil et al. 2006a b). These high velocities fa¬ 
vor a scenario where the MCs are in their first infall 
towar d the MW ( [Besla e t al.|2007 Boylan-Kolchin et al. 


2011 


iBusha et aT7 


2011 h In this scenario, Besla et al. 


12012) found that the Magellanic Stream may have been 


formed by LMC tides on the SMC before the system 
was accreted by the MW. Diaz & Bekki (2012) studied 
several orbits consistent with the recent proper motions. 
They assumed an orbital history where the MCs became 
a strongly interacting pair just recently. In this picture, 
their first close encounter provided enough tidal forces 
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to disrupt the SMC’s disk and create the Magellanic 
Stream. In contrast to the Magellanic Stream, t he Mag- 


ellanic Bridge has a young stellar population (Irwin et al. 
19851, besides being composed of gas and dust (|Hind- 

-„ I i .,.11 i nn a l llonn.u 1/ r — J — 


man et al. |1963[ |Sofiie||T994] |Muller et al.||2Q0i] | (Jordon 
et al. 1200*5 ). The study of gas-to-dust ratio of the Magel¬ 


lanic Bridge confirms the hypothesis that the Magellanic 
Bridg e is fo rmed of materi al that was stripped from the 
SMC (Gordon et al. 20091. Moreover it was found that 
at the Magellanic Bridge the gas-to-dust ratio is about 
12 times higher than in the Galaxy and about two times 
higher than the value expected by its metalicity, indicat¬ 
ing that some of the dust in that region must have been 
destroyed. 

It is not well known whether the MCs were formed as a 
bina ry system or were dyna mically coupled about 3 Gyr 


ago (Bekki & Chiba 


2005). The kinematic history re¬ 


const nrcnoiToTTEieTMCssuggests that their last closest 
approach occurred about 0.2 Gyr ago, when they came 
within 2-7 kpc distance of each other, a likely progenitor 
event for the construction of the Bridge. Another pos¬ 
sible result of this interaction is the stretchi ng in north¬ 


easte rn and Wing sections toward the LMC (Westerlund 
19911. 

The last interaction was also likely resp onsi ble f or the 


bimodality of the SMC along its depth (Mathewson & 
Ford 1984), a feature that is corroborated by the ve- 


locity com ponents in HI data, which are 40 km s 1 in 


difference (|Mathewson 1984 Stanimirovic et al. 2004|). 
Distances to Gepheids show that the high velocity com¬ 
ponent is located further in distance, while the low veloc¬ 
ity component is closer. Both components appear to be 
expanding at around 15 km s -1 . The low velocity com¬ 
ponent itself is further subdivided into two other groups, 
one extending from the center to the northeast (NE) at 
a location of 50 kpc, the other, extending from the Bar 
to the southwest (SW) at a location of 60 kpc. The Bar 
is the most extended region (45-90 kpc) and therefore 


the HI profiles in th at direction are very complex (Math- 
ewson fe Ford|[l984). The stud y of the HI component a t 


the SMC shows that it rotates (Stanimirovic et al.|2004), 
indicating that the SMC contains a disk component. 

Distances to Cepheids also show that the two gas com¬ 
ponents are about 10 kpc apart and have a depth of 
17 kpc. A double component, relating to the distances, 
is also observed in Cepheids located in the Bridge, con¬ 
sistent with the two gaseous components. Besides the 
Bridge’s bifurca tion, it is also, on avera ge, closer to us 


than the SMC (Mathewson et al. 


1986). Recently, red 


clump (RC) and RR Lyrae stars (RRL S) were used to 
estimate the SMC’s line of sight depth. |Subramanian fe| 
Subramaniam (20121 used these stars to study the three- 
dimensional structure of SMC and found that SMC has 
an average depth of ~ 14 kpc. They concluded also that 
the NE region is closer to us and an elo ngation along the 
NE-SW direction is seen. |Nidever et al. (2013|) used RC 
stars to study the eastern and western regions of SMC 
and found that the first has a larger depth (~ 23 kpc), 
while the second has a much shallower depth (~ 10 kpc). 
Some of their eastern fields showed a bimodality, with one 
component located at ~ 67 kpc and the other at ~ 55 kpc 
away from us, leading them to conclude that the closer 
stellar component was stripped from SMC during the 


past interaction. The study of the stellar distribution of 
SMC shows that it has a sp heroid al or sli ghtly ellipsoidal 
structure (Subramanian & Subramaniam 2012). A stel¬ 
lar disk structure is also observed; however, there is a 
discre pancy between the gaseous and stellar d i sk param¬ 
eters (j Grocncwcgcn 2000[|Dobbie et al. ||20 1 U |Subrama- 
nian & Subramaniam 2(jl5). 

The overall picture shows that SMC possesses a stellar 
and a gaseous disk component, as well as a spheroidal 
stellar component. The eastern region and Magellanic 
Bridge are closer in distance to us and show a bimodal 
structure. These regions are also the ones that were prob¬ 
ably most affected by the past interaction, which explains 
why they may be approaching LMC, which is located at 
~ 50 kpc from us. The line of sight depth estimations 
show a relative high dispersion from region to region, 
with the eastern side being deeper than the western side. 

1.1. The SMC’s Magnetic Field 

The magnetic field of SMC was previously studied us¬ 
ing several techniques - optical interstellar polarization, 
Faraday rotation, and synchrotron emission. The inter¬ 
stellar polarization is caused by the alignment of the dust 
grains’ angular momentum with a local magnetic field. 
The unpolarized light emitted by the stars become polar¬ 
ized due to the dichroic extinction by dust grains. There¬ 
fore, optical polarization vectors trace the magnetic field 
projected in the plane of the sky, and a polarization map 
can be interpreted as a magnetic field map in that direc¬ 
tion. The phenomenon known as Faraday rotation hap¬ 
pens, when radiation travels through an ionized medium 
with a local magnetic field. The light plane of polar¬ 
ization rotates, if the local magnetic field has a parallel 
component to the direction of propagation of the light. 
The incident light is decomposed in to two parts circu¬ 
larly polarized with opposite rotations. The synchrotron 
emission is produced when charged particles are accel¬ 
erated by a magnetic field. The radiation emitted by 
these particles is polarized, because the acceleration is 
n ot isotro pic. 

Mao et al. (|2008[) u sed opt ical polarization data from 


Mathewson & Ford (1970a) and obtained an ordered 


magnetic held projected in the plane of the sky of 
B s k y = (1.6 ± 0.4) ii G and a turbulent component 
of 2 [iG. Magalhaes et al. (2009) a nalyze d optical polar¬ 
ization data from Magalhaes et al. (1990) and estimated 
an ordered sky-projected m agnetic field of 1. 7 fiG and 
a turbulent field of 3.5 /uG. Mao et al. (2008) also esti¬ 
mated the line of sight ordered magnetic held, through 
Faraday rotation, using data from the Australia Tele¬ 
scope Compact Array (ATCA). They obtained Blos = 
(0.19 ± 0.06) fj,G, with the field pointing away from us. 
In addition, there are several studies based on the ra- 
dio continuum synch r otron emission 


Loiseau et al. 1987 


(Haynes et al.|1986 
Haynes et al. 1990p. Using syn¬ 


chrotron emission and assuming energy equipartition be¬ 


tween the cosmic rays and magnetic fields, Loiseau et al. 
(1987) estimated the SMC’s total magnetic held to be 
about 5 gt,G. 


Pho t oelectric polarization dat a from Magalhaes et al. 
dl99 0|); IMagalhaes et al. | (|2009|) |S chmidt| (1970111976), 
and |Mathewson &; Ford| (|1970a|b|) suggested that SMC 
has an ordered magnetic held aligned with the Magel¬ 
lanic Bridge. Polarized synchrotron emission, despite its 
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weak nature i n the SMC, also confirmed such an align¬ 
ment (Hayn es et al. 1990|). Besides the alignment with 
the Bridge, |Magaihacs et al.| (|2009|) observed an align¬ 
ment of the ordered magnetic field with the SMC’s Bar. 
The magnetic field may have played an important role in 
shaping the Bridge, since the interstellar magnetic field 
of galaxies can be strong enough to influ ence the dynam¬ 
ics of their gas (Zweibel & Heiles 1997). Therefore, this 
work aims to study the structure of the magnetic field 
using CCD data, which is more precise than photoelec¬ 
tric. We concentrated our study on the NE and Wing 
regions of the SMC, since these regions may have been 
more affected by the last interaction with LMC. We also 
have some data of the Magellanic Bridge, from the part 
closer to the SMC, a structure probably formed due the 
interaction of the SMC-LMC system. 

The paper is structured as follows. In Section [2j we 
describe the observational data, the reduction process, 
and how we built the catalog. The estimates for the 
foreground polarization are presented in Section [3] We 
analyze the polarization patterns for each field observed 
in Section HI The magnetic field geometry is discussed 
in Section 15] The alignment of the ordered magnetic 
field with respect to the Magellanic Bridge is discussed 
in Section [6] We estimate the magnetic field strength in 
Section]?] Finally, in Section[8] we summarize our results 
and the main conclusions. 


2. OBSERVATIONS AND DATA REDUCTION 
2.1. Observational Data 

The observational data were taken at Cerro Tololo 
Inter-American Observatory (CTIO) over the course of 
five nights in 1992 November 13-17, applying the op¬ 
tical polarization technique and using the 1.5 m tele¬ 
scope. The instrumental setup consisted of a half-wave 
retarder plate followed by a calcite, V filter, and the de¬ 
tector. This setup gives us, for each star, two images 
with orthogonal polarizations. At least two sets of im¬ 
ages, taken with the half-wave retarder plate positioned 
at different angles, are necessary for each object/f ield , to 
obtain the linear Stokes’ parameters Q and U (Serkowski 
1974). Some standard stars polarization were measured 
using 16 images (see Table [If. Further details about the 
instru mental technique can "be found in|Magalhac s et al.| 
(1996). A TeklK-1 lk x lk CCD, with a plate scale of 
0".434/pixeL read noise of 8.0602 e^ rms, and gain of 
9.005 e“/ADU was used. 

During the observational run, bias and flat-field im¬ 
ages were taken to correct instrumental noise. One non¬ 
polarized standard star was observed, in order to check 
whether there was any instrumental polarization or not. 
Also, three polarized standard stars were observed, in 
order to get the conversion for the polarization angles 
(6»n into the equatorial system. Twenty-eight fields of 
8 x8 arcmin size, situated at the NE and Wing sections 
of the SMC as well as the Magellanic Bridge, were ob¬ 
served. We decided to concentrate our data efforts on 
these regions for the reasons mentioned previously, relat¬ 
ing to the last int eract ion with the LMC, about 0.2 Gyr 
ago ( Westerlund||1991 ). 


6 The polarization angle 9 is most common refereed as the po¬ 
sition angle (PA). Throughout this paper we use the notation 9 
instead of PA. 


The date of the observation, number of positions the 
retarder plate was rotated, and integration time for the 
standard stars are presented in Table [l] The coordinates 
and date of the observation for the 28 fields are presented 
in Table [2j For the fields, the half-wave retarder plate 
was rotated four times and the integration time was 300 s 
per frame. Figure [I] shows the positions of these fields. 


Table 1 

Observations of standard stars 


HD 

Date 

No. a 

IT b 

(s) 

Standard 

9540 

1992 Nov 13 

16 

20 

Non-polarized 

283812 

1992 Nov 13 

16 

8 

Polarized 

23512 

1992 Nov 14 

16 

3 

Polarized 

23512 

1992 Nov 15 

4 

3 

Polarized 

23512 

1992 Nov 16 

4 

3 

Polarized 

23512 

1992 Nov 17 

4 

4 

Polarized 

298383 

1992 Nov 17 

4 

6 

Polarized 


a Number of images. 
b Integration time in each image. 


Table 2 

Observations of SMC fields 


SMC a 

R.A. b 

(h:m:s) 

Decl. b 

Date 


01 

01:00:50.2 

-71:51:55.2 

1992 

Nov 

13 

02 

01:00:49.8 

-72:06:56.2 

1992 

Nov 

13 

03 

01:00:49.4 

-72:21:55.2 

1992 

Nov 

13 

04 

01:00:49.0 

-72:36:56.2 

1992 

Nov 

13 

05 

01:03:48.9 

-71:51:57.0 

1992 

Nov 

13 

06 

01:03:48.5 

-72:06:58.0 

1992 

Nov 

13 

07 

01:05:16.9 

-72:37:27.8 

1992 

Nov 

14 

08 

01:06:46.7 

-72:21:59.8 

1992 

Nov 

14 

09 

01:08:15.6 

-72:36:59.7 

1992 

Nov 

14 

10 

01:08:15.1 

-72:52:00.7 

1992 

Nov 

14 

11 

01:09:44.2 

-72:59:31.7 

1992 

Nov 

14 

12 

01:12:41.7 

-73:29:33.7 

1992 

Nov 

14 

13 

01:14:11.9 

-73:07:03.7 

1992 

Nov 

14 

14 

01:15:40.6 

-73:22:04.7 

1992 

Nov 

14 

15 

01:17:10.2 

-73:14:35.8 

1992 

Nov 

14 

16 

01:21:39.4 

-72:44:39.1 

1992 

Nov 

17 

17 

01:21:38.2 

-73:14:40.1 

1992 

Nov 

15 

18 

01:24:35.8 

-73:37:11.4 

1992 

Nov 

16 

19 

01:30:32.4 

-73:52:17.3 

1992 

Nov 

15 

20 

01:42:26.8 

-73:52:26.9 

1992 

Nov 

15 

21 

01:45:23.3 

-74:30:00.8 

1992 

Nov 

15 

22 

01:48:23.7 

-74:00:03.7 

1992 

Nov 

15 

23 

01:51:22.8 

-73:52:36.7 

1992 

Nov 

16 

24 

02:01:46.6 

-74:15:17.8 

1992 

Nov 

15 

25 

02:00:15.8 

-74:37:45.2 

1992 

Nov 

17 

26 

02:06:13.1 

-74:37:52.9 

1992 

Nov 

15 

27 

02:09:12.8 

-74:22:54.3 

1992 

Nov 

15 

28 

01:54:37.9 

-74:30:20.1 

1992 

Nov 

15 


a Field’s label. 
b Coordinates in J2000. 


2.2. Reduction Process 

The data were reduced using the software IRABU more 
specifically the NOAO and PCCD packages. This last 
package w as developed b y the Polarimetric Group from 
IAG/USP ( Pereyra|2000 ) to deal with polarimetric data. 


7 IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Observed SMC fields overlapped with a Spitzer/MIPS image at 160 //m ((Gordon et ah 2011). The square size represents the 
real FOV. The fields with no polarization trend are labeled as P0; the fields with one polarization trend are labeled as PI; the fields with 
two polarization trends are labeled as P2; the fields that were filtered by magnitudes are labeled as PFM; and the fields that had stars 
excluded are labeled as PSE. The classification for the fields is introduced at Section 111 


We eliminated instrumental noise applying overscan, 
bias, and flat-field corrections. Once the images were 
processed, the Stokes Q and U parameters were obtained 
performing aperture photometry of ordinary and extraor¬ 
dinary images for each stellar object in the field. This 
was done in an automated manner, firstly looking for all 
objects with a stellar profile, then performing the pho¬ 
tometry for each of these, finally using the set of im¬ 
age pairs the Stokes parameters Q and U were obtained. 
M ore deta ils about the polarimetric tasks can be found 


Pereyra (2000). 

A polarization of P 


= (0.035 ± 0.027)% was obtained 
for the non-polari zed star HD9540, a value within lcr of 
the published one (Heiles|2000). The instrumental polar¬ 
ization was therefore found to be negligible. For the po¬ 
larized stars, we found polarization int ensities, w i th val ¬ 
ues within 3cr of on es publis hed by Hsu fe Breger (1982), 
Bastien et al. (1988), and HPOL (2011). One exception 
to this was HL)298383, that lay within ~ 3.4cr. The dis¬ 
crepancy for this measurement is d ue to the very small 
error of the published value by Tapia (1988). The good 


agreement indicates that the instrument measured the 
polarizations precisely. Lastly, using the polarization an¬ 
gles from the polarized standard stars, we estimated the 
Ad conversion for the equatorial system. The standard 
deviation for A 9 is 1°.8, indicating a good determination 
of the conversion factor. 

After the reduction process, we carefully analyzed the 
stars that had suspiciously high polarizations. Stars 
too bright, too faint, or with badly centered coordinates 
tended to have such high polarizations. Since the mea¬ 
surements for these stars are not reliable, we discarded 
them. 


2.3. Polarimetric Catalog 

The polarimetric catalog contains 7207 stars and was 
constructed using stars with P/crp > 3. This new catalog 
is a huge improvement in the number of stars com pared 
to previous polarimetric catalog s for the SMC (|Mathew- 
son & Ford 1970a Schmidt 1976). It consists of positions, 
polarization intensity and its associated error, polariza¬ 
tion angle, and V magnitude and the corresponding er¬ 
ror. A short version of the polarimetric catalog, for the 
observed and intrinsic polarizations, can be seen in the 
Appendix [A] 

In order to obtain the position of the stars in right 
ascension (R.A.) and declination (deck), we made use 
of Digital Sky Survey (DSS) images. These images had 
the same center and size as our fields. Then, using a 
reference star, we performed a transformation from pixel 
to astronomical coordinates. The average errors for the 
positions are ap.A. = 0.25 s and <Td ec i. = 1.0". About 
12% of the stars have repeated coordinates, due to the 
close proximity of a bright and faint star in the images. 
This occurs in the astrometry procedure when the faint 
stars are re-centered at the bright star positions, due to 
software limitations. We nevertheless decided to keep 
those stars in the catalog since the exact positions are 
not important for the further magnetic field study. 

The V magnitude (my) was obtained through sum¬ 
ming up the flux of the image pair using one frame. In 
order to convert these magnitudes into the UBV system, 
we made use of published catalogs. Tabic [3] presents the 
references for these catalogs and the number of stars used 
to perform the calibration for each field. For some fields, 
we used more than one star to estimate the Amy cali¬ 
bration parameter. For others, using the same star we 
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looked for its my in different catalogs. The average er¬ 
ror for my is a mv = 0.13 mag, which was obtained by 
averaging the individual errors. These errors take into 
account the photon noise error and the dispersion in the 
determination of Amy. 

To evaluate the completeness of the catalog, we con¬ 
structed a histogram with bin sizes of 1 mag. The com¬ 
pleteness or limit magnitude was determined by the his¬ 
togram’s peak. Afterwards, we evaluated the standard 
deviation in the bin of the limiting magnitude. A new 
histogram with bin sizes equal to the standard deviation 
( a mv =0.28 mag) was constructed, which can be seen 
in Figure [2] Finally, the new limiting magnitude was de¬ 
fined as the new peak of the histogram. Using this pro¬ 
cedure, the catalog’s completeness is my = 18.00 mag. 
We suspect that most of the stars with my < 14 mag 
belong to the Galaxy. In Section [3] we further discuss 
this assumption. 



10 12 14 16 18 20 22 

n v (mag) 


Figure 2. V magnitude histogram of the catalog stars. The ver¬ 
tical line at my = 14.2 mag separates the SMC stars from the 
foreground candidates. It is also indicated the limit magnitude 
{my lim = 18.00 =b 0.28 mag) of the catalog. 


3. FOREGROUND POLARIZATION 

The light from the SMC’s member stars has to travel 
through the Galaxy’s ISM. Therefore, part of the polar¬ 
ization we measure is due to the MW’s dust. The SMC’s 
interstellar polarization is obtained by subtracting the 
foreground polarization from the Galaxy, which is esti¬ 
mated in this section. 


Schmidt (1976) divided the SMC into five regions (see 
Figure 2 in this paper) and estimated the foreground po¬ 
larization for these regions using foreground stars. |Ro-| 
drigues et al.| (|1997|) re-estimated the foreground polar¬ 
ization tor each of these regions and obtained more pre¬ 
cise values. In our analysis, we used the same regions 
from Schmidt (1976) and added a sixth region between 
the SMC’s body and the Magellanic Bridge, which we 
call region IV-V. One of the motivations to estimate the 
foreground polarization from our data is to be able to 
subtract the Galactic foreground in that new region. 

Our polarimetric data is spread along the NE and Wing 
regions of SMC and the Magellanic Bridge. Therefore it 
does not include objects from Schmidt (1976)’s region I 


and that is why we do not have an estimate for the Galac¬ 
tic foreground there. Region II corresponds to the NE 
part of the SMC. Regions III and IV represent the SMC’s 
Wing, with the former located nearer to the SMC’s Bar. 
Region IV-V is located between regions IV and V, the 
latter being located furthest relative to the SMC and at 
the Magellanic Bridge. 

Our catalog has both SMC and Galaxy stars. Thus, 
one way to estimate the foreground polarization is to de¬ 
termine which stars are located in the Galaxy and use 
their polarizations as a measure for the foreground. Be¬ 
low we estimate a magnitude threshold for SMC stars. 

The distance modulus relation is 


TOy — My = 5 log d — 5 + Ay, 


( 1 ) 


where my is the apparent magnitude in the V band; My 
is the absolute magnitude in the V band; d is the distance 
in pc; and Ay is the MW foreground extinction in the V 
band. 


The SMC’s distance is d = (63 ± 1) kpc (Cioni et al. 
20001, which was obtained usin g the tip of the r ed gy 
ant branch method. Recently, Scowcroft et al. (20151 
obtained d = (62 ± 0.3) kpc, using Spitzer observations 
of classical Cepheids, the new measurement agrees with 
the distance we assumed. In order to estimate the aver¬ 
age foreground extinction, one can use the average f ore¬ 
g roun d color excess, E(B — V) ~ 0.05 mag (Bessell 
19911, together with the av erage extincti on law for the 
MW, Ay ~ 3.2 E(B — V) ([Seaton 1979]), thus we get 
Ay ~ 0.2 mag. Hence, the distance modulus to the SMC 
is my — My ~ 19.2 mag. 

The brightest giant and super-giant stars have 
My < — 5 mag, which translates into my < 

14.2 mag when considering t he previo usly mentioned dis¬ 
tance modulus for the SMC. Massey (2002)’s catalog fo¬ 
cused on a precise magnitude calibration for bright stars 
from the SMC, to complement other catalogs focused on 
faint stars. This catalog has 84995 stars with my < 
18 mag, from which just 3.5% have my < 14.2 mag. 
Therefore, considering that all stars with my < 14.2 mag 
are Galactic objects, there is a possible inclusion of only 
a small number of SMC members. 


In addition, a l l the stars used by |Schmidt| (fl976| and 
Rodrigues et al.| (|1997|) to estimate the foreground polar¬ 
ization have my W 14.2 mag. Considering this fact and 
the argument above, it is a reasonable approximation to 
characterize the stars from our catalog with magnitudes 
below this value as foreground star candidates. Using 
this limit we would be excluding from the foreground 
sample just the Galactic late-type dwarfs. 

A second filter was applied to select objects representa¬ 
tive of foreground polarization. We excluded all objects 
with polarizations higher than 1% and differing from the 
average, of the stars with my < 14.2 mag, by a factor 
larger than 2cr. This procedure may remove stars with 
intrinsic polarization. 

The steps to separate the foreground stars and esti¬ 
mate the foreground polarization were the following: 


1. Take the stars with my < 14.2 mag and P < 1%; 


2. Obtain the Stokes parameters Q and U and its as¬ 
sociated error for the selected stars; 
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Table 3 

Photometric calibration 


SMC a 

No. of Stars b 

Reference 0 





01-19 

4 

Massey11|2002| 

20-21 

1 

Ardeberg 6z, Maurice (1977) ,|Isserstedt (1978), and Lasker et al. (2007J 

22 

1 

isserstedt ( i9/'8) and Ardeberg 

(1980) 

23 

2 

.Lasker et ai. 

2u07 

) and 

banduJ 

elk (1989) 

24 

1 

.Lasker et ai. 

2007 

i and 

Pene et al. (1991) 

25 

1 

Lasker et al. 

2UU7 




26-27 

3 

Demers 77 Irwin 1x99111 



28 

2 

kasker et al. 

(2O0Y 





a SMC field labels. 
b Number of stars used. 

c Reference for the catalogs used to make the my calibration. 


3. Divide the stars into the Schmidt (1976)’s regions 
and region IV-V; 


4. Obtain the uncertainty weighted average Stokes pa¬ 
rameters for each region and its standard deviation; 


5. Exclude the stars with Stokes parameters deviating 
more than 2cr from the average; 


6. Using the remaining stars, the Stokes parameters 
are re-computed using an uncertainty-weighted av¬ 
erage. We estimated the error of the Galactic 
polarization in each region as the mean standard 
deviation, because the error from the uncertainty 
weighted average is usually underestimated; 

7. Finally, from the Stokes parameters, we obtain the 
polarization, its angle, and the error for each re¬ 
gion. 


Table [4] s hows our estim ate f or the foreground alon g 
side that of|Schmidt| (|1976|) and |Rodrigues et al.| (|1997|) . 
It can be seen that our estimates are in good agreement 
with the two previously mentioned works. We obtained 
smaller errors due to the higher precision of our CCD 
measurements compared with previous photoelectric po- 
larinretry. In addition, the ratio of foreground polariza¬ 
tion to foreground extinction is smaller than 3 % ma g^ 1 
in all regions, as expected for the Galaxy (W hittet|l 992). 
These facts provide a reasonable degree ot confidence" in 
our estimates. We used our values to subtract the fore¬ 
ground polarization from the stars of our catalog. The 
foreground subtraction was done via the Stokes parame¬ 
ters. 

Figure [3] shows histograms for the ratio between the 
intrinsic and observed polarizations ( Pint/Pobs ) and the 
difference between them (0j nt — 9 0 b s ). The foreground 
subtraction can modify the polarization intensities up to 
50% from the observed value. Likewise, the polarization 
angles can shift up to 20°. There is no large difference 
between the different foreground corrections. 


4. POLARIZATION PATTERNS IN THE SMC 
FIELDS 

With the aim of analyzing the magnetic field geometry, 
we studied the polarization preferential direction (PD) 
for each one of the 28 fields. We used stars with intrinsic 
polarization P/crp > 3, subtracting the foreground us¬ 
ing the estimate obtained in this work. The number of 
objects in each field is shown in Table [5j 


In order to obtain the PDs from 9 , we performed Gaus¬ 
sian fits to the 9 histograms. To find the polarization 
intensity (P) for the PDs, we took the field stars with 
9 within 3 a of the mean 9 obtained by the Gaussian 
fits. The polarization intensity for the PDs was esti¬ 
mated from the median of the polarization for these stars. 
We used the median instead of the average because it is 
a better indication of the central trend for asymmetri¬ 
cal distributions. This procedure was not done for fields 
with a random distribution in 9. 

Our fields showed different distributions for the 9 his¬ 
tograms, therefore we classified the fields into five differ¬ 
ent groups: P0, PI, P2, PFM, and PSE. Table [5] presents 
the classification for each field. The next Sections explain 
what they are and how we dealt with them. In order to 
test in a more robust way our field’s classification, we 
performed the P-test on our sample. The results were 
consistent with our previous classisfication. More details 
about the P-test procedure can be found in the Appendix 
|B| Plots of the polarization map, 9 vs my, 9 histogram, 
and pola rization intensity histogram are shown in Fig¬ 
ures 4.1-4.28 (online material). 


4.1. Fields with No PD (P0) 

The fields SMC 05, 14, 19, and 25 did not show any 
PD, therefore we did not fit any Gaussian to them. Fig¬ 


ure |4.1 shows the example of SMC 14: all the stars in 
the field were considered to construct the 9 and P his¬ 
tograms. This field shows a completely random distribu¬ 
tion. Hence, we are not able to define a ny P D for the 
polarization. In the case of SMC05 (Figure [Ot), there is 
possibly a PD; however, the number of stars is not high 
enough to characterize it. These fields are represented 
by gray squares in Figure [I] 

4.2. Fields with One PD (PI) 

The 9 histogram for fields SMC 03, 06, 09, 11, 15, 18, 
21, and 28 showed one PD. Consequently one Gaussian 
was fitted for each of these. So me fi elds show a clear PD, 
an example is SMC 28 (F igure fo) ;). In others, SMC09 
for instance (Figure |4.4( :), it is almost possible to see a 
second PD; however, the number of stars is once again 
not sufficient for the Gaussian fitting. Finally, some fields 
dis playe d a large dispersion, an example is SMC03 (Fig¬ 
ure 


4.5:), but not large enough to qualify the distribu¬ 


tion’s nature as random. The red squares in Figure |T] 
represent these fields. 

4.3. Fields with Two PDs (P2) 
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Table 4 

Foreground polarization for each region 


Region 3, 

Pfor,B b 

(%) 

0for,s C 

(deg) 

Pfor,r b 

(%) 

A C 

u for,r 

(deg) 

-^for.lg b 

(%) 

0for,lg C 

(deg) 

Nt i d 

JV for,lg 

Schmidt (19761 

Rodrigues et al. (1997 1 

This work 


i 

0.37 ± 0.15 

111 

0.47 ± 0.09 

113.6 

— 

— 

— 

ii 

0.27 ± 0.15 

123 

0.30 ±0.08 

124.1 

0.316 ±0.016 

111.6 

75 

hi 

0.06 ± 0.09 

139 

0.17 ± 0.05 

145.0 

0.124 ±0.024 

136.8 

12 

IV 

0.14 ±0.12 

125 

0.22 ±0.05 

124.2 

0.162 ± 0.020 

125.6 

16 

IV-V 

- 

- 

- 

- 

0.138 ±0.048 

102.7 

15 

V 

0.16 ±0.12 

93 

0.27 ± 0.05 

95.6 

0.301 ± 0.028 

91.1 

20 


3 Regions’s label. 

b Foreground polarization intensity (Pf or x)- 
c Foreground polarization angle (0f or ,x)- 

d Number of stars used for the foreground estimation by this work. 




Figure 3. Histograms of Pint/P^h * and (6,nt — On h*)- The red histogram is obtained after 
orange histogram is obtained after Rodrigues et al. (1997} foreground subtraction, and the 
subtraction from this work. The dashed lines represent no foreground. 


|Schmidt| ([1976} foreground subtraction, the 
blue histogram is obtained after foreground 
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Figure 4.1. Plots for SMC14, an example of a field with no PD (labeled as PO). Panel (a) shows a polarization map overlapped with a 
Spitzer /MIPS image at 160 pm ( 1 Gordon et al.|[~2011| ); panel (b) shows a plot of 0 versus magnitude; panel (c) shows a 6 histogram; and 
panel (d) shows a polarization intensity histogram. 


SMC05 



AR (J2000) m v (mag) 6 (deg) P (%) 


Figure 4.2. Plots for SMC05, an example of a field with no PD (labeled as P0). Panel (a) shows a polarization map overlapped with a 
Spitzer/MIPS image at 160 pm (I Gordonet al. |[20111 ); panel (b) shows a plot of 9 versus magnitude; panel (c) shows a 8 histogram; and 
panel (d) shows a polarization intensity histogram. 


Two PDs for fields SMC 01, 02, 04, 07, 08, 10, 22, 
and 27 were observed in the 9 histogram. Similarly, two 


Gaussians were fitted to reproduce the two PDs. Figure 
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Figure 4.3. Plots for SMC28, a n example of a field with one PD (labeled as PI). Panel (a) shows a polarization map overlapped with 
a Spitzer /MIPS image at 160 //m i|Gordon el al. 2011); panel (b) shows a plot of 0 versus magnitude; panel (c) shows a 8 histogram, the 
green line represents the Gaussian fit; and panel (dj snows a polarization intensity histogram, the vertical line on this plot represents the 
median value. 



Figure 4.4. Plots for SMC09, a n example of a fie ld with one PD (labeled as PI). Panel (a) shows a polarization map overlapped with 
a Spitzer/MIPS image at 160 /mi ijGorcion et al. 201 lj); panel (b) shows a plot of 0 versus magnitude; panel (c) shows a 8 histogram, the 
green line represents the Gaussian fit; and panel (d) snows a polarization intensity histogram, the vertical line on this plot represents the 
median value. 



Figure 4.5. Plots for SMC03, a n example of a fie ld with one PD (labeled as PI). Panel (a) shows a polarization map overlapped with 
a Spitzer/MIPS image at 160 /mi l|Gordon et al. 2011); panel (b) shows a plot of 0 versus magnitude; panel (c) shows a 8 histogram, the 
green line represents the Gaussian lit; and panel (dj snows a polarization intensity histogram, the vertical line on this plot represents the 
median value. 


|4.6| : shows the 9 histogram for SMC01, which, along with 
other fields, a possible third PD can be discerned. Gaus- 
sians were just fitted to those PDs that are significant 
enough relative to the random background and distin¬ 
guishable from the neighboring PD. In the other hand, 
some fi elds display two clear PDs, an example is SMC10 
(Figure 4.7 z). The magenta squares in Figure [l] represent 
these fields. 


4.4. Fields Filtered by Magnitudes (PFM) 

The 9 histogram for fields SMC 12, 13, 20, and 23 
are the result of a PD and a random background being 
superposed. In these fields, the peak position was well 
fitted, but the standard deviation was not. To improve 
the fits, we made plots of 9 vs my, to check whether it 
is possible t o se parate the random background from the 
PD. Figure 4.8 shows, for SMC12, the histograms for 9 
and P, as well as the plot for 9 vs my. For this specific 


case, the random background is due to stars with my > 
17.5 mag. Hence, we selected the stars with magnitudes 
smaller than this value to perform the Gaussian fit. This 
cut in magnitude can be justified, taking into account 
that magnitudes are distance indicators, therefore the 
objects forming the random distribution may constitute 
a further population. In order to determine the optimum 
value of my, trying to separate the random background, 
we used plots of the dispersion of 9 in function of my. 
More details about this procedure can be seen in the 
Appendix [C] These fields are represented by steel blue 
squares in Figure [T] 


4.5. Fields with Stars Excluded (PSE) 

For the remaining fields, SMC 16, 17, 24, and 26, a 
simple Gaussian fit does not converge if all objects are 
considered. In the case of SMC 24, there are two PDs in 
the 9 vs my plot (Figure 4.9 d). One of them is com- 
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Figure 4.6. Plots for SMC01, an example of a f ield with two PDs (labeled as P2). Panel (a) shows a polarization map overlapped with 
a Spitzer /MIPS image at 160 firn i|Gordon el al. 2011); panel (b) shows a plot of 8 versus magnitude; panel (c) shows a 8 histogram, the 
green and gold lines represent the Gaussian fits tor tne first and second PDs, respectively; and panel (d) shows a polarization intensity 
histogram, the vertical lines on this plot represents the median value for the first (solid line) and second PD (solid-dotted line). 




Figure 4.7. Plots for SMC10, an example of a fi eld with two PDs (labeled as P2). Panel (a) shows a polarization map overlapped with 
a Spitzer/MIPS image at 160 fim |]Gordon et al. 2011); panel (b) shows a plot of 8 versus magnitude; panel (c) shows a 8 histogram, the 
green and gold lines represent the Gaussian nts for tne first and second PDs, respectively; and panel (d) shows a polarization intensity 
histogram, the vertical lines on this plot represents the median value for the first (solid line) and second PD (solid-dotted line). 



Figure 4.8. Plots for SMC12, an example of a field for which we had to filter the stars by magnit udes to fit the polar ization PD (labeled 
as PFM). Panel (a) shows a polarization map overlapped with a Spitzer/MIPS image at 160 (im l|Gordon et al.|2011(l; panel (b) shows a 
plot of 8 versus magnitude, the vertical line demarcates the magnitude cut; panel (c) shows a 6 histogram, the green line represents the 
Gaussian fit; and panel (d) shows a polarization intensity histogram, the vertical line on this plot represents the median value. 


posed of few stars, such that two Gaussian fits were not 
possible. Hence, we decided to exclude these stars. For 
the other three fields, we could not observe any differ¬ 
ent behavior al ong th e magnitude range, an example is 
SMC16 (Figure |4.10|}). In order to be able to perform 
the Gaussian fits, we arbitrarily excluded the stars that 
were preventing the fit convergence. The spring green 
squares in Figure [l] represent these fields. 


5. MAGNETIC FIELD GEOMETRY 

In this Section we discuss the magnetic field geometry 
of the SMC based on polarization maps. Initially, a pre¬ 
liminary check was made by constructing a 9 histogram 
in order to check the general behavior of the data. We 
used the data that was foreground-corrected with the 
foreground estimate from this work (see Section [3]) and 
considered objects with intrinsic polarization: P/ap > 3. 
Figure [5] shows the 0 histogram. 


The histogram in Figure [5] reveals three major trends 
in the sample. The most prominent is featured at 
(65° ± 10°) (trend I). The latter two trends were identi¬ 
fied at (115°±10°) (trend II) and (150°±10°) (trend III). 
The errors for each trend were estimated as the bin size 
considered for the histogram. To assess whether these 
trends are either separated in distance or if any kind of 
segregation is present, density plots (Q vs. I/, 6 vs. P, 
9 vs. my, and P vs. my) were constructed, as can be 
seen in Figure [6] These trends were further reinforced 
in the plots of Q vs. U , 9 vs. P, and 9 vs. my. We 
do not observe any strong segregation; however, trend 
II is concentrated at brighter stars (in the range ~ 17.2 
to ~ 18.25 mag) and smaller polarization intensities (in 
the range ~ 0.25% to ~ 0.75%). This indicates that 
trend II could be associated with the part of the bimodal 
structure that lies closer to us in distance. Trend I en¬ 
compasses a wide range in polarization (from ~ 0.2% to 
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SMC24 



AR (J2000) m v (mag) 9 (deg) P (%) 


Figure 4.9. Plots for SMC24, an example of a field for which we had to exclude some stars to fit the polar ization PD (labeled as PSE). 
Panel (a) shows a polarization map overlapped with a Spitzer/MIPS image at 160 pm (Gordon et al. |201l|; panel (b) shows a plot of 6 
versus magnitude, the stars in the light pink area were not considered for the Gaussian tit; panel (c) shows a 6 histogram, the green line 
represents the Gaussian fit; and panel (d) shows a polarization intensity histogram, the vertical line on this plot represents the median 
value. 
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Figure 4.10. Plots for SMC16, an example of a field for which we had to exclude some stars to fi t the polar ization PD (labeled as PSE). 
Panel (a) shows a polarization map overlapped with a Spitzer/MIPS image at 160 pm (Gordon et al.||2011|; panel (b) shows a plot of 0 
versus magnitude, the stars in the light pink area were not considered for the Gaussian fit; panel (cj shows a 6 histogram, the green and 
gold lines represent the Gaussian fits for the first and second PDs, respectively; and panel (d) shows a polarization intensity histogram, the 
vertical lines on this plot represents the median value for the first (solid line) and second PD (solid-dotted line). 



6 (deg) 


Figure 5. 0 histogram of the data corrected by foreground from 
this work and with P/crp > 3. The vertical lines demarcate the 
trends. 

~ 0.9%) and magnitude (from ~ 17.1 to ~ 18.9 mag), 
while trend III is shifted in the vertical direction in both 
polarization (from ~ 0.4% to ~ 0.8%) and magnitude 
(from ~ 17.5 to ~ 19 mag), but concentrated in a smaller 
range of polarization intensities. In the P vs. my plot 
no segregation is observed, just the normal behavior of 
polarization increasing with magnitude. 

In order to check the impact of the foreground sub¬ 
traction on the data, we made density plots of the same 
quantities as before, but using the SMC observed polar¬ 
ization (Figure [7]). For the observed data, just trend II 


is visible, which demonstrates how much the foreground 
subtraction changes the geometry observed, raising two 
other trends that were masked by the foreground contri¬ 
bution. 

Based on the discussion above, we separated the PDs 
obtained in the previous section in to three groups: 

1. 40° < 6> pd < 90° (trend I); 

2. 90° < 6> pd < 132°.5 (trend II); 

3. 132°.5 < 6> pd < 185° (trend III). 

The upper and lower limits for each group were chosen 
to be the position of its trend plus/minus half distance to 
the neighbor trend. Trend I has no neighbor in the left 
side, therefore its lower limit was defined as the position 
of trend I minus half distance between trend I and II. 
Trend III has no neighbor in the right side, therefore its 
upper limit was defined as a value that includes all the 
PDs lying in its right side. 

Each group contains one of the trends, thus all the PDs 
were classified as belonging to a trend. Figure [8] shows 
the maps obtained for each group and the foreground 
map as determined by this work. The foreground map 
(Figure^) shows that the Galactic foreground is highly 
aligned with the Bridge direction at a PA of 115°.4. Fig¬ 
ure shows that the fields in the NE region present an 
ordered magnetic field roughly aligned with the Bar at a 
PA of ~ 45°. Figure [8]: also shows an alignment with the 
Bridge direction; however, it is the trend containing the 
least number of vectors. Figure [8}i shows vectors sepa¬ 
rated by about 90° from the ones in Figure [8jo. Overall 
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e (deg) P (%) 

Figure 6. Density plots of the SMC intrinsic polarization. Panel (a) shows a plot of U vs. Q, panel (b) of P vs. 6, panel (c) of my vs. 
9, and panel (d) of my vs. P. The color bars indicate the number of objects. 


the polarization maps do not show a very demarcated 
structure, but rather a quite complex geometry. 

5.1. Magnetic Field Origin 

The origin of t he SMC ’s large scale magnetic field was 
discussed by Mao et al. (20081. This study concluded 
that a cosmic-ray driven dynamo can explain the exis¬ 
tence of a large scale magnetic field at the SMC in terms 
of time scale arguments. Nonetheless, it had difficul¬ 
ties in explaining the geometry observed, because of the 
unidirectionality of the magnetic field and no change of 
sign for the Faraday rotation measures (RM). Consid¬ 
ering that the cosmic-ray driven dynamo is the mech¬ 
anism that generated the large scale magnetic field at 
the SMC, initially this field had mostly an azimuthal 
configuration, which is the axisymmetric dynamo mode 
m = 0. Tidal interactions lead to the excitation of the 
bisymmetric mode m = 1 , when the axisymmetric mode 
is already at work (Moss 1995). The tidal interactions 
also induce nonaxisymmetric velocities in the interact¬ 
ing galaxies disks and may lead to the damping of the 
m = 0 mode, leaving just the bisymmetric magnetic field 
( Vogler fe Schmitt||2001 ). 


Mao et al. (2008) explained that in the scenario above, 
the tact that the SMC’s magnetic field is predominantly 
lying on its disk, can be explained by the fact that the az¬ 
imuthal magnetic field produced by a dynamo mostly lies 
in the disk of the galaxy. The inclusion of tidal forces be¬ 
tween the SMC and LMC would explain the alignment 
with the Magellanic Bridge. Finally, if the SMC pos¬ 
sesses a bisymmetric magnetic field, we would observe a 
periodic double change of the RM sign with respect to the 
azimuthal angle. If the magnetic field is represented by 
a superposition of m = 0 and m = 1 modes, even more 
sign changes would be expected. Considering that the 
magnetic field lines do close, but the locations with field 
lines pointing toward us are outside the SMC’s body, this 
would explain why it is observed just negative RMs. The 
regions that should display positive RMs may have a low 
emission measure of ionized gas, therefore RM is zero in 
this locations, since to observe non-zero RM the average 
electron density in the line of sight should be non-zero. 
Our d ata do n ot exc lude the physical explanation given 
by Mao et al. (2008). On the contrary, the unidirection¬ 
ality tor the magnetic field is no longer a problem. 

Trend I is widely spread in magnitudes and polariza- 
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0 (deg) P (%) 

Figure 7. Density plots of the SMC observed polarization. Panel (a) shows a plot of U vs. Q , panel (b) of P vs. 9, panel (c) of my vs. 
9, and panel (d) of my vs. P. The color bars indicate the number of objects. 


tions, which may in dicate its correlation with a global 
pattern of the SMC. Stanimirovic et al. (2004) obtained 
that the PA for the major kinematic axis of the SMC is 
around 50°, which is 15° in difference from trend I. This 
direction is also roughl y the PA of the SMC’s Bar, po¬ 
sitioned at ~ 45° (van den Bergh 2007). This suggests 
that in the Bar region the magnetic - field may be cou¬ 
pled to the gas and therefore the field lines are roughly 
parallel to the Bar direction due to the flux freezing con¬ 
dition; however, we can not explain why the field lines 
in the SMC’s Wing and in the Magellanic Bridge display 
also an alignment with this direction. Our understand¬ 
ing is that the initial m = 0 dynamo mode may have 
been damped by the nonaxisymmetric velocities excited 
by the tidal interactions. Therefore we do not observe the 
typical behavior of an azimuthal field in the polarization 
vectors. Nonetheless, the bisymmetric field was left and 
possibly higher order dynamo components. The current 
geometry for the SMC’s magnetic field is probably the 
product of an active interaction of the SMC-LMC-MW 
system summed to the influence of star formation, su¬ 
pernova explosions, and other processes that can inject 


energy to th e ISM, controlling its dynamics. |Bur khart 
et al. (2010) suggests that the HI in the SMC is super- 
hic 


super- 

Aifvemc, which also explains the rather disordered con¬ 


figuration observed for the large scale magnetic field. 

As mentioned before, our data are concentrated at the 
NE and Wing sections of the SMC and at a part of the 
Magellanic Bridge. In these regions a bimodality in the 


ewson et al. 1986 Nidever et al. 


velocity~cbmponcnts are observed 


distance is known to exist from Cepheid distances (|Math- 

1 


2013[). Similar ly, two 


ewson 


studies (Math- 


1 1981 iStanimirovic et al.|2004[ ). Figure [6] indicates 
rend II could be related to the component located 


that trend II could be related to the component 
closer to us and that trend III could be linked to the 
most distant component. The difference in magnitudes 
between trends II and III is about 0.5 mag. Translating 
into relative distances: dni = 1.26c?//, neglecting inter¬ 
nal extinction. Hence, if trend II belongs to the closest 
component, located a t 55 kpc as observed from Cepheids 
in the eastern region (Nidever et al. 2013]), trend III may 
be located at 69 kpc. Considering the photometric er¬ 
ror of our catalog of 0.13 mag, the distance for trend III 
obtained by this rough estimate is compatible with the 
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Figure 8. Polarization maps. Panel (a) shows a map of the foreground polarization from this work, panel (b) shows a map of the PDs in 
the range [40°.0 — 90°.0], panel (c) shows a map of the PDs in the range [90°.0 — 132°.5], and p anel i d ) shows a map of the PDs in the 
range [132°.5 — 185°.0]. The vectors are overlapped with a S'pitzer/MIPS image at 160 //m ([Gordon et. al.|20l ! ). The Magellanic Bridge 
stars in the lower left region of the maps and extends up to the LMC location. 


distance obtained by Cepheids for the furt hest compo¬ 
nent, which is 67 kpc (Nidever et al. 2013). The tidal 
interactions between SMC and LMC are likely to ex¬ 
plain the stretching of the magnetic field lines toward 
the Magellanic Bridge direction. Nonetheless, this ef¬ 
fect was important just in the component closer to us 
in distance, which is also closer to the LMC. The cre¬ 
ation of bridges and the magnetic field alignment with 
respect to the brid ge was alread y observed in numeri¬ 
cal studies. Kotarba et al. (2011) simulated the interac¬ 
tion of three disk galaxies up to the point where they all 
merge. Their simulation shows that the magnetic field 
of the interacting galaxies strongly changes with time ac¬ 
cording to their interaction. Nonetheless, the comparison 
of the SMC-LMC-MW system with their results is not 
straightforward, because the properties of the galaxies 
are different. 

The coincidence between the directions of trend II and 
the Galactic foreground is not easy to explain. The mag¬ 
netic field at the Galactic halo is not expected to be high 
and indeed the polarization measurements are rather low 
in that region (P/ or < 0.5%). A possible speculation is 
that the MW halo also feels the tidal forces by the SMC 
and LMC, therefore its magnetic field is also stre tched 
in the same direction. The simulation by Kotarba et al. 


(2011) shows that when the galaxies are about 50 kpc 
apart their magnetic fields align in the outskirts of the 
approaching galaxies (Figure 15 in their paper). The 


masses, magnetic fields strengths, and 3D distribution 
of the SMC-LMC-MW system is not alike their system. 
Nonetheless, we speculate that the coincidence of trend 
II with the Galactic foreground can be explained by the 
system interaction. 

Yet, a question can be raised regarding the genuinity 
of trend II: is it real or a vestige of a bad foreground 
removal? This question is difficult to address since the 
mentioned direction is that of the Galactic foreground. 
If the foreground was underestimated, the trend could be 
a remnant of the foreground itself; however, the median 
polarizations for trend II (0.46% < P/j < 1.2%) 

range two to eight times above the estimated foreground 
(0.06% < Pf or < 0.47%). Therefore a foreground 
underestimation can be justifiably ruled out. Another 
possibility is that faint stars from the MW are included in 
the catalog, causing the foreground trend to persist. To 
fully rule this out, the distance to those stars or another 
distance indicator such as the E(B — V ) color excess 
is necessary to separate SMC and MW members. The 
upcoming GAIA mission will prove to be a good tool 
for SMC-MW member separation due to the parallaxes 
that will be measured in the MW halo. The accuracy 
expected for the fainter stars (my = 20 mag) is to be 
as good as 1 mas, about 10% for stars at 10 kpc. The 
stars that form trend II have magnitudes from 17.2 to 
18.3 mag, thus the GAIA accuracy may be good enough 
to define whether these stars belong to the MW’s halo 















14 


Lobo Gomes, Magalhaes, Pereyra & Rodrigues 


or the SMC. 

This work brings a new understanding of the SMC’s 
magnetic field. Neverth eless, the SOUTH POL project 
(Magalhaes et al.||2012 1 will measure the polarization of 
objects in the whole southern sky (south of —15° ini¬ 
tially), which will increase even more the polarization 
sample toward the SMC. These new data will help to 
get a more complete picture of the SMC magnetic field 
structure, since it will measure polarizations in regions 
not covered by this work. Usi ng the photometric cata¬ 
log of SMC members of Massey (2002), to get the num¬ 
ber of stars per magnitude range, we expect that polar¬ 
izations will be measured for around 7500 bright stars 
(my < 15 mag) with accuracy up to 0.1% and around 
68,000 stars (my < 17 mag) with accuracy up to 0.3%. 
Naturally, fainter stars will de detected, but with accura¬ 
cies that may or may not be appropriate for ISM studies; 
for instance, stars with my ~ 18 mag will be measured 
with about 1% of accuracy. Moreover, SOUTH POL will 
measure additional foreground objects toward the SMC, 
which should have GAIA distance estimates. 

6. ALIGNMENT BETWEEN THE SMC MAGNETIC 
FIELD AND THE MAGELLANIC BRIDGE 

To further address the question regarding the align¬ 
ment between the SMC polarization and Magellanic 
Bridge direction, we used cumulative frequency distri¬ 
bution (CFD) analysis. Firstly, we calculated the angle 
between the SMC and LMC centers: 9 qm = 115°. 4, 
which is in the sa me direct ion as the Magellanic 
Bridge. Similarly to [Schmidt (19761, we defined 9 m 
as the angles between #om and the polarization angles 
of the stars. The following coordinates were used 
for the centers of the SMC and LMC respectively: 
R.A.(J2000) = 00 /t 52 m 38.0 s , decl.(J2000) = -72°48'oi" 
and R.A.(J2000) = 05 ,l 23 m 34.6 s , 

decl.(J2000) = —69° 45'22". 

Considering stars with P/ap >3 and my > 14.2 mag, 
CFDs were constructed. Table [5] shows the number of 
stars used per held. Figure [9] shows the results for six sets 
of data: (a) using all the stars from the catalog, (b)-(f) 
using stars from region II to region V. For each set of data 
four CFDs were evaluated: using the observed polariza¬ 
tion and the foreground-corrected polarization consider¬ 
ing the three aforementioned estimates. The straight line 
at 45° represents a random distribution. If a distribution 
is above this line, 9m is concentrated at smaller values, 
which indicates an alignment with the SMC-LMC direc¬ 
tion. A distribution below the straight line indicates the 
magnetic held being perpendicular to the SMC-LMC di¬ 
rection. 

For all cases, barring region III, the 9m distribution 
for the observed polarization is above the straight line, 
indicating an alignment with the SMC-LMC direction. 
The foreground-corrected CFDs, considering all the dif¬ 
ferent estimates, lie below the observed one. This correc¬ 
tion brings the distributions closer to the straight line. 
This demonstrates the importance of a proper foreground 
subtraction to understand a possible SMC magnetic held 
alignment with the SMC-LMC direction. This is particu¬ 
larly true because the foreground polarization is approx¬ 
imately aligned with the SMC-LMC direction (Figure 
go- For all the cases, except region IV-V, the distribu¬ 
tions roughly follow the straight line until a noticeable 


deviation at around 40°. The location of the deviation 
corresponds to the direction of trends I and III, as pre¬ 
viously mentioned, at 65° and 150° (Oqm +40 — 90 and 
9qm + 40). As briefly touched upon, region III does not 
follow the above explanation. In this case, the distri¬ 
butions, including the observed polarization, lie slightly 
below the straight line. Similarly to the other regions a 
noticeable deviation up from the straight line can be seen 
at around 30°. Region IV-V shows a similar behavior, 
but only for the foreground-corrected distribution. 

In order to check whether we made an appropriate 
choice for the definition of 9qm, we repeated the anal¬ 
ysis this time adopting 9qm as the angle between the 
star’s position and the LMC’s center. Overall there is no 
qualitative difference in the results. 

To quantify by how much our CFDs a re sim ilar to_the 
unifor m one, we used the approach from |Rodrigues et al.| 
(2009). This involved doing a Kuiper test, which is a 


variant of the Kolmogorov-Smirnov test and more ap¬ 
propriate for cyclic quantities such as 9m- Most of the 
distributions had probabilities much smaller than 1 % of 
being uniform , exce pt for: (i) region IV using the Ro- 
drigues et al.| (|1997D foreground estim ate, the probabil¬ 
ity is of 34%; (h) region V using the |Rodrigues et al. 
(1997) and this work foreground estimates, for which the 
probabilities are 8% and 7%, respectively. The probabil¬ 
ities are not very high, therefore none of the cases can be 
classified as a uniform distribution. 

Previous studies observed an overall alignment of the 
SMC’s ordered magnetic field and the Magellanic Bridge 
direction. Our sample shows such overall alignment only 
for the observed data. Once the foreground is removed 
a complex geometry for the magnetic field arises. The 
CFDs analysis displays neither a uniform distribution 
nor a strong alignment with the Bridge direction. The 
SMC’s magnetic field, despite being complex, is not to¬ 
tally random. The new strong alignments correspond to 
trend I and III, while trend II is not predominant in the 
CFDs study. 

7. MAGNETIC FIELD STRENGTH 

Our data can also be used to estimate the magnetic 
field streng th. For this purpose, we used the |Chan- 


drasekhar fc Fermi (1953) method modified by Falceta- 


Gongaives et al. "j 21)1)8). Since the interstellar polariza¬ 


tion is due to an alignment of the dust grains’ angular 
momentum with a local magnetic held, one can expect 
that for a larger magnetic held, the dispersion of the po¬ 
larization angles will be smaller. 

Besides the polarization angles’ dispersions (68), which 
were already estimated using the Gaussian hts, we also 
need to know the gas velocity dispersion in the line of 
sight direction (SVlos) and the gas mass density (p), in 
order to be able to apply this method. These quantities 
were estimated using m aps of HI velocity d i spersion and 
HI column density from Stanimirovic et al. (1999 2004). 


7.1. HI Velocity Dispersion 

In our analysis we used a HI velocity dispersion map 
that was derived using combined images from the ATCA 
and Parkes. The second-momentum analysis was used to 
obtain the velocity dispersion map, w hich can be seen in 
Figure 4 of Stanimirovic et al. (2004). 
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B u (deg) 

Figure 9. Cumulative frequency distribution for 6m , the angle between the polarization angles and the SMC- LMC direction. The plot on 
the upper left, panel (a), shows the CFD for all stars; on the upper middle, panel (b), for region II; on the upper right, panel (c), for region 
III; on the bottom left, panel (d), for region IV; on the bottom middle, panel (e), for region IV-V; and on the botto m right, panel (f), 
for region V. 0 m,o6s corresponds to the observed polarization, 6m, fs corresponds to polarization corrected by |Schmidt| (1976|l foreground, 
6M,fr corresponds to polarization corrected by|Rodrigues et al.|(1997| foreground and 6M,fq corresponds to polarization corrected by 
foreground from this work. The dot-dashed lines correspond to a uniform distribution. 


Since this map includes just the SMC’s body, SVlos 
was locally estimated just for those fields coinciding with 
the map, i.e., SMC01-19. The estimate for each field 
was obtained by averaging the line of sight velocities 
in the area of the SVlos map coinciding exactly with 
each of our 8x8 arcrnin SMC fields. For fields lo¬ 
cated at the Magellanic Bridge, we cons i dered the av¬ 
erage value for this region from Bruns et al. (2005): 
SVlos = 20kms _1 . 


7.2. HI Mass Density 

The same data set of Section 17.11 was used. The HI 
velocity profiles were integrated to obtain the HI co lumn 
density map, which is shown in Figure 3 of|Stanimirovic 
eFahl|l999|. 

"We obtained the local estimates as previously de¬ 
scribed for fields SMC01 -19. For the remaining fields, we 
used the average v alue of N(HI) = 5 x 10 20 atoms cm -2 
(Brims et al .|2005 ) for the Bridge. Before averaging over 
the squares to calculate N(HI), the aforementioned quan¬ 
tity was corrected by a factor to account for the hydrogen 
auto-absorption along the line of sight and defined by the 


following equation: 

'1 + 0.667 {log N(HI) - 21.4); for logN(HI) > 21.4, 


/ = 


1; for logN(HI) < 21.4. 


( 2 ) 

This correction is required for regions with column den¬ 
sities higher than 2.5 x 10 21 atoms cm -2 (Stanimirovicj 
etH1[l999t . 

In order to convert from HI column density [ N(HI )] 
to HI number d ensi ty (nr/), the SMC’s d epth had to be 
estimated. Subramanian & Subramaniam (2012) used 
the dispersions in the magnitude-color diagrams of RC 
stars together with distance estimates of RRLS to de¬ 
termine the SMC’s depth along the line of sight. They 
obtained an average depth of (14 ± 6) kpc (error ob¬ 
tained by private communication with the first author). 
As mentioned in Section [TJ the line of sight depth at 
SMC varies from region to region, therefore adopting an 
average value leads to some uncertainty in the determi¬ 
nation of uh- For the fields located at the Wing and 
Bridge, this average is a good esti mate. For the fields a t 
the NE, where the depth is higher ( Nidever et al.||2013 ), 
we may be underestimating this value, which leads to 
an overestimation of njj- Nevertheless, the mass den- 
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sity was obtained by applying p = 7 mnun ■ Considering 
the SMC’s abundances the e quivalent molecular weight 
is 7 = 1.22 ( Mao et aL|2008 ). 


7.3. Magnetic Field Strength Estimates 


Following Falceta-Gongalves et al. (2008) method, we 
firstly determined the turbulent magnetic field strength 
(SB), assuming equipartition between the turbulent mag¬ 
netic field energy and turbulent gas kinetic energy (Equa- 
tion|3|. Lastly, we determined the ordered magnetic field 
strength on the plane of the sky (B s k y ) through Equation 


SI 


\pWlos - ^ B 2 , 


(3) 

(4) 


Table [5] shows the obtained magnetic field intensities 
for each field, as well as the values for the HI velocity dis¬ 
persion at the line of sight and HI number density. For 
the fields with more than one PD, different estimates 
were calculated using the different polarization angles’ 
dispersion, which are also shown in Table [5j When we 
obtain more than one PD for the polarization angle dis¬ 
persion, we know that they must be located at different 
distances or direction and for the first case the values of 
SVlos and nn would be more appropriate for the most 
distant PD. It is not easy to quantify the uncertainty for 
our estimates, but we want to point out that the errors 
might be quite high. Nevertheless, the usage of the inte¬ 
grated values lead to an overestimation of uh and SVlos > 
therefore overestimating SB and underestimating B s k y . 

The fields SMC 03, 06, and 18 have negative estimates 
for B s ky The fact that we may be underestimating B s k y 
can be an explanation for the negative values; nonethe¬ 
less, these fields are also the ones with larger SO, all of 
which are higher than 7 t/4 . For dispersions of that order, 
the magnetic field may have just a turbulent component 
or a high in clinati on wit h respe ct to the plane of the sky 
(Falceta-Gongalves et al.] |2008 ). We know from Faraday 
rotation that there is an ordered line of sight magnetic 
field of B los = (0.19 ± 0.06) p G ( |Mao et al.|[2008] ), 
despite this small value there is the possibility that in 
these regions Blos dominates over B s k y . It should also 
be noted that all the values of B s k y (with exception of the 
negative ones), and certainly the average B s k y value, are 
significantly higher than the average Blos- This shows 
that the magnetic field of the SMC is, in general, mostly 
in the plane of the sky. This lends additional value to 
the study of the B s k y structure in the SMC. 

The Bsk y dispersion is quite high as can be seen in 
Figure |T0j In many cases, the values can be up to 25 
times above the average. This high dispersion is not so 
surprising, because our estimates may contain not just 
measurements of the diffuse interstellar magnetic field, 
but also estimates for local structures (e.g., shells, clus¬ 
ters, molecular clouds, HII regions). Optical polarimetric 
observations already showed that some ISM structures 
are magnetized, possessing magnetic fi elds of up to tens 
of pG, for i nstance, IRAS Vela Shell (Pereyra fe M ag- 
alhaes|2007 ) and NGC 2100 ( Wisniewski et al.||2007 1, so 
that some of the high values of B s k y could be associated 
to such structures. To verify this hypothesis we queried 


SimbacJ^] in the regions of our fields and looked for ISM 
structures. Table [6] presents a list of objects per field. 
Some of the measurements may be related to these struc¬ 
tures; nonetheless, we can only guarantee that there is 
spatial projected correlation between the structures and 
the polarization vectors. With respect to the fields that 
possess magnetic fields of tens of pG, SMC04, for in¬ 
stance, possesses a shell with a size of 5x5 arcmin, this 
structure might be responsible for our estimated value 
for H s ky = 23.6 pG. The fields with a large SO, SMC 06 
and 14, for instance, are the fields with a largest number 
of structures. It is possible to find direct associations of 
the polarization angles with the geometry of these struc¬ 
tures, but for doing this a detailed study of each source 
would be necessary, which is beyond the scope of this 
work. 


20 


15 


5 


0 


Figure 10. Histogram of the ordered sky-projected magnetic 
field. The vertical line at 0.947 fiG represents the uncertainty 
weighted average (excluding the negative values). 


< B sky > 
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A turbulent magnetic field value of 

SB = (1.465 ± 0.069) pG was obtained from 

the uncertainty weighted average, when considering 
all fields. The same analysis, including both observed 
components, led to the computation of the ordered 
magnetic field projected on the plane of the sky: 
Bsky = (0.947 ± 0.079) pG. The negative values were 
not used to evaluate the average. Our estimate for SB 
is about 30% and 60% smaller than the value obtained 


by Mao et al. (2008) and Magalhaes et al. (2009) 
respectively. The trend of smaller values is also obs erved 
in the results for B s k y , which is 40% smalle r than |Mao 


et al.| (|2008) and 50% reduced relative to Maga) 


haes 


(2009). The high difference in these estimates arise 


mam 


et al. 


y due to t he way nn was evalua ted. Both Mao| 
(20081 and Magalhaes et al. (2009) used a constant 


value for the HI number density (nn = 0.1 atoms cm -3 ). 
We can see in Table [5] that the fields SMC20-28 have 
nn one order of magnitude smaller than the value 
quoted in the previously mentioned papers, therefore 
our estimates for the magnetic field in these regions 
should naturally lead to smaller values, also reducing 
the average. The observation of a weak large scale 
magnetic field agrees with the expectation that the HI 


http://simbad.u-strasbg.fr/simbad/ 
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Table 5 

Local ISM parameters for our fields 


SMC a 

SVLOS h 

n H 

r c 

se d 

SE 

> e 

D f 

J-'sky 

06 

P h 

Label 1 

Nod 


(km s" 

- 1 ) 

(atoms cm 3 ) 

(deg) 

(M g) 

(mG) 

(deg) 

(%) 



01 

12.643 ± 

0.078 

0.102 ± 

0.043 

21.28 ±0.29 
22.9 ±1.2 

2.05 ± 

0.44 

3.21 ± 0.69 
2.80 ±0.66 

120.90 ±0.28 
46.3 ± 1.1 

0.90 

0.85 

P2 

491 

02 

14.763 ± 

0.091 

0.092 ± 

0.039 

40.93 ± 0.90 
19.5 ± 2.2 

2.27 ± 

0.48 

0.35 ±0.11 
4.1 ± 1.2 

58.39 ± 0.73 
150.0 ± 1.6 

0.62 

0.66 

P2 

492 

03 

24.228 ± 

0.075 

0.144 ± 

0.061 

46.08 ±0.61 

4.66 ± 

0.99 

—0.17 ± 0.10 

102.16 ±0.55 

1.0 

PI 

471 

04 

25.723 ± 

0.072 

0.183 ± 

0.078 

39.7 ±1.9 

10.8 ±2.2 

5.6 ± 

1.2 

1.14 ±0.51 
23.6 ± 7.9 

169.2 ± 1.6 
65.5 ±2.0 

0.87 

1.0 

P2 

170 

05 

13.799 ± 

0.093 

0.174 ± 

0.074 

- 

2.92 ± 

0.62 

- 

- 

- 

P0 

47 

06 

19.690 ± 

0.085 

0.164 ± 

0.070 

71.9 ±3.2 

4.04 ± 

0.86 

-2.72 ±0.63 

100.6 ±2.1 

1.2 

PI 

247 

07 

27.824 ± 

0.045 

0.165 ± 

0.070 

43.0 ±2.2 
25.28 ± 0.72 

5.7 ± 

1.2 

0.42 ± 0.48 
6.4 ± 1.4 

84.7 ± 1.6 
178.85 ± 0.73 

1.0 

1.0 

P2 

640 

08 

22.09 ± 

0.13 

0.104 ± 

0.044 

37.03 ±0.51 
21.76 ±0.88 

3.61 ± 

0.77 

1.18 ±0.26 
5.4 ± 1.2 

73.69 ± 0.38 
154.79 ±0.66 

0.87 

0.93 

P2 

735 

09 

26.332 ± 

0.066 

0.154 ± 

0.066 

18.84 ±0.16 

5.2 ± 

1.1 

10.1 ± 2.2 

67.51 ± 0.16 

0.98 

PI 

605 

10 

23.764 ± 

0.063 

0.154 ± 

0.066 

26.45 ± 0.36 
22.67 ± 0.72 

4.7 ± 

1.0 

4.8 ± 1.0 

6.6 ± 1.5 

142.54 ±0.31 
54.97 ± 0.60 

0.94 

1.0 

P2 

560 

11 

22.443 =fc 

0.071 

0.160 ± 

0.068 

17.17 ±0.17 

4.55 ± 

0.97 

10.2 ±2.2 

158.29 ±0.17 

1.4 

PI 

453 

12 

22.200 ± 

0.049 

0.083 ± 

0.035 

14.55 ± 0.68 

3.24 ± 

0.69 

9.3 ±2.1 

147.88 ± 0.68 

0.74 

PFM 

205 

13 

19.89 ± 

0.14 

0.163 ± 

0.069 

24.1 ± 2.6 

4.07 ± 

0.86 

5.0 ± 1.5 

126.2 ± 2.5 

0.46 

PFM 

122 

14 

21.08 ± 

0.14 

0.137 ± 

0.058 

- 

3.95 ± 

0.84 

- 

- 

- 

P0 

165 

15 

23.31 ± 

0.11 

0.134 ± 

0.057 

27.2 ± 1.1 

4.32 ± 

0.92 

4.08 ± 0.96 

144.5 ± 1.0 

1.1 

PI 

139 

16 

20.77 ± 

0.12 

0.056 ± 

0.024 

25.6 ±4.6 
22.8 ±7.0 

2.48 ± 

0.53 

2.7 ± 1.2 

3.4 ±2.1 

134.5 ± 4.8 
66.8 ±6.4 

0.89 

0.72 

PSE 

87 

17 

20.602 ± 

0.069 

0.097 ± 

0.041 

30.0 ±3.2 
15.0 ±1.6 

3.26 ± 

0.69 

2.38 ±0.89 
8.9 ±2.3 

118.9 ±2.6 
54.0 ± 1.6 

0.87 

0.92 

PSE 

136 

18 

20.842 ± 

0.084 

0.059 ± 

0.025 

59.5 ± 3.7 

2.57 ± 

0.55 

-1.06 ±0.32 

88.5 ±2.9 

1.1 

PI 

131 

19 

18.96 ± 

0.14 

0.046 ± 

0.019 

- 

2.05 ± 

0.44 

- 

- 

- 

P0 

40 

20 

20 


0.0115 ± 

0.0049 

20.6 ±3.2 

1.09 ± 

0.23 

1.80 ±0.62 

75.0 ±3.2 

0.69 

PFM 

63 

21 

20 


0.0115 ± 

0.0049 

7.51 ± 0.42 

1.09 ± 

0.23 

7.1 ± 1.6 

159.57 ±0.49 

1.7 

PI 

57 

22 

20 


0.0115 ± 

0.0049 

35.2 ± 5.4 

21 ± 15 

1.09 ± 

0.23 

0.45 ± 0.32 
1.7 ±2.2 

157.8 ±4.3 

75 ± 14 

0.96 

0.89 

P2 

60 

23 

20 


0.0115 ± 

0.0049 

8.9 ± 1.3 

1.09 ± 

0.23 

5.8 ± 1.6 

158.8 ± 1.2 

0.75 

PFM 

63 

24 

20 


0.0115 ± 

0.0049 

27.9 ±2.4 

1.09 ± 

0.23 

0.96 ±0.29 

150.1 ± 2.4 

1.3 

PSE 

63 

25 

20 


0.0115 ± 

0.0049 

- 

1.09 ± 

0.23 

- 

- 

- 

P0 

40 

26 

20 


0.0115 ± 

0.0049 

31.3 ±3.2 

1.09 ± 

0.23 

0.70 ±0.27 

79.4 ±2.9 

0.92 

PSE 

60 

27 

20 


0.0115 ± 

0.0049 

22.7 ±2.9 
7.8 ± 1.2 

1.09 ± 

0.23 

1.51 ± 0.49 
6.8 ± 1.9 

69.4 ±2.6 
116.2 ± 1.0 

0.69 

0.67 

P2 

72 

28 

20 


0.0115 ± 

0.0049 

13.90 ±0.77 

1.09 ± 

0.23 

3.30 ± 0.74 

69.87 ± 0.77 

0.82 

PI 

72 


a Field’s label. 
b HI velocity dispersion. 
c HI number density. 

d Polarization angle dispersion, for fields with two PDs the first row shows the first PD and the second row the second PD. 
e Turbulent magnetic field. 

f Ordered sky-projected magnetic field, for fields with two PDs the first row shows the first PD and the second row the 
second PD. 

g Trend for the polarization angle, for fields with two PDs the first row shows the first PD and the second row the second 
PD. 

h Median polarization, for fields with two PDs the first row shows the first PD and the second row the second PD. 

1 Field’s classification. 

J Number of stars with P/<rp > 3. 


in the SMC is super-Alfvenic ( Burkhart et al.|[2010 ). 

The average turbulent component is higher than the 
average ordered component, a common result for all 
types of galaxies. The B s k y /5B = 0.65 ratio for the J3MC 
is closer to the MW, ranging from 0.6 to 1.0 (Beck|2001), 
than the typical va lues for other irregular dwarfs, ~ 0.2 
(Chyzy et al. 2011). The production of magnetic fields 
in irregular dwarf galaxies is most likely not maintained 
by a large-scale dynamo pr ocess, due to the small rates 
observed for B s k y /SB (Chyzy et al. 2011). In the case 
of the SMC, this process can be the mechanism that cre¬ 
ates the large-scale magnetic field observed in the SMC, 
as discussed in Section [5] 


8. SUMMARY AND CONCLUSIONS 

This work used optical polarimetric data from CTIO, 
aiming to study the magnetic field of the SMC, an irreg¬ 


ular galaxy and satellite of the MW. One of the biggest 
peculiarities of the SMC is that its ISM is particularly 
different from that of the Galaxy (e.g., high gas-to-dust 
ratio and submm excess emission), most likely due to the 
large difference in metalicity. The data reduction led to 
a catalog with 7207 stars, with well determined polariza¬ 
tions (P/ap > 3). This new catalog is a great improve¬ 
ment compared to previous catalogs for the SMC. 

Our analysis showed that caution is necessary when 
subtracting the foreground Galactic polarization, be¬ 
cause this correction strongly changes the geometry ob¬ 
served for the magnetic field. We present a new estimate 
for the foreground Galactic polarization using the stars 
from our catalog, which has smaller errors compared to 
previous ones. 

This catalog was used to study the magnetic field on 














18 


Lobo Gomes, Magalhaes, Pereyra & Rodrigues 


Table 6 

ISM structures in our fields 


SMC a 

Object Name 

R.A. b 

(h:m:s) 

Decl. b 

Object Type 

01 

DEM S 107 

01:00:08.7 

-71:48:04 

HII region 


[MOH2010] NE-lc 

01:00:10.0 

-71:48:40 

molecular cloud 


DEM S 113 

01:01:30.3 

-71:47:45 

HII region 


[SSH97] 290 

01:01:12 

-71:52.9 

shell 


LHA 115-N 72 

01:01:32.5 

-71:50:43 

HII region 

02 

[SSH97] 279 

01:00:26 

-72:05.1 

shell 


[SSH97] 285 

01:00:46 

-72:09.6 

shell 

03 

[SSH97] 288 

01:01:09 

-72:24.6 

shell 

04 

[JD2002] 21 

01:01:33.6 

-72:34:53 

planetary nebula 


[SSH97] 273 

01:00:05 

-72:39.7 

shell 

05 

DEM S 121 

01:03:03 

-71:53.8 

HII region 


[SSH97] 319 

01:04:40 

-71:53.8 

shell 

06 

[MOH2010] NE-3g 

01:03:10.0 

-72:03:50 

molecular cloud 


DEM S 119 

01:03:01.9 

-72:05:17 

HII region 


LHA 115-N 76B 

01:03:08.0 

-72:06:24 

HII region 


LHA 115-N 76A 

01:03:48.9 

-72:03:52 

HII region 


[SSH97] 309 

01:03:03 

-72:03.2 

shell 


SNR B0101-72.6 

01:03:17 

-72:09.7 

supernova remnant 


MCELS S-148 

01:03:48.6 

-72:03:56 

HII region 


MCELS S-147 

01:03:25.0 

-72:03:45 

HII region 


MCELS S-144 

01:03:01.2 

-72:05:41 

HII region 

07 

[SSH97] 329 

01:05:18 

-72:34.4 

shell 


[SSH97] 330 

01:05:23 

-72:38.9 

shell 


[SSH97] 321 

01:04:42 

-72:37.5 

shell 

08 

DEM S 134 

01:06:48.5 

-72:24:32 

HII region 


[SSH97] 344 

01:06:46 

-72:25.0 

shell 


[SSH97] 355 

01:07:30 

-72:21.8 

shell 

09 

[SSH97] 367 

01:08:43 

-72:38.1 

shell 


[SSH97] 358 

01:07:33 

-72:39.9 

shell 

10 

[SSH97] 354 

01:07:25 

-72:54.9 

shell 


DEM S 133 

01:07:34.7 

-72:51:20 

HII region 

11 

[SSH97] 376 

01:09:29 

-72:57.4 

shell 


[SSH97] 386 

01:10:31 

-73:03.1 

shell 


[SSH97] 373 

01:09:20 

-73:01.9 

shell 


[SSH97] 387 

01:10:33 

-73:01.6 

shell 

12 

[BLR2008] SMC N83 4 

01:12:05.8 

-73:31:01 

molecular cloud 


[SSH97] 398 

01:12:23 

-73:28.0 

shell 


[BLR2008] SMC N83 2 

01:12:41.3 

-73:32:14 

molecular cloud 

13 

[SSH97] 408 

01:13:30 

-73:03.6 

shell 


[SSH97] 407 

01:13:29 

-73:03.6 

shell 

14 

2MASX J01144713-7320137 

01:14:47.132 

-73:20:13.80 

planetary nebula 


DEM S 152 

01:14:54.1 

-73:19:45 

HII region 


NAME SMC B0113-7334 

01:14:44.9 

-73:20:06 

HII region 


DEM S 157 

01:16:20 

-73:20.2 

HII region 


[SSH97] 423 

01:15:29 

-73:23.8 

shell 


MCELS S-193 

01:14:55.7 

-73:20:10 

HII region 


MCELS S-195 

01:15:04.7 

-73:19:10 

HII region 


MCELS S-191 

01:14:41.7 

-73:18:06 

HII region 

15 

[SSH97] 440 

01:17:27 

-73:12.5 

shell 


DEM S 159 

01:16:58 

-73:12.1 

HII region 


DEM S 155 

01:17.1 

-73:14 

HII region 

17 

LHA 115-N 87 

01:21:10.69 

-73:14:34.8 

planetary nebula 


[SSH97] 472 

01:21:40 

-73:18.0 

shell 

19 

[MSZ2003] 28 

01:30:44 

-73:49:42 

shell 


[MSZ2003] 32 

01:31:43 

-73:52:24 

shell 

20 

[MSZ2003] 49 

01:42:36 

-73:49:54 

shell 


[MSZ2003] 48 

01:41:35 

-73:55:12 

shell 

21 

[MSZ2003] 52 

01:45:23 

-74:30:06 

shell 


[MSZ2003] 55 

01:46:10 

-74:28:24 

shell 

26 

[MSZ2003] 88 

02:06:37 

-74:36:48 

shell 

27 

[MSZ2003] 95 

02:09:22 

-74:24:48 

shell 


the SMC. After foregroun^ d removal, J ?fir^'e trends at the following polarization angles were observed: (65° ± 10°), 
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(115° ± 10°), and (150° ± 10°). For the first trend, the 
polarization vectors in the NE region are roughly aligned 
with the Bar direction, which is at a PA of ~ 45°, rein¬ 


forcing what was observed by Magalhaes et al. (2009). 
In the case of the second trend, the polarization angle 
is aligned with the Bridge direction, which is at a PA of 
115°.4, and possess the same direction as the Galactic 



A possible explanation tor the magnetic held alignment 

with the Bar direction is that the magnetic field is cou¬ 
pled to the gas, therefore the field lines are parallel to the 
Bar direction due to the flux freezing condition. The sec¬ 
ond trend may be due to tidal stretching of the magnetic 
field lines in the direction of the Magellanic Bridge. The 
coincidence of the alignment of the second trend with 
the Galactic foreground may be due to the Galactic halo 
also feeling the tides from the MCs in the region close 
to the MCs, therefore the MW’s halo magnetic field gets 
also stretched in the same direction. The third trend 
does not display any particular feature. Regardless the 
trends, the magnetic field structure seems to be rather 
complex in the SMC. 

The polarization and magnitude distributions of the 
115° and 150° trends suggest that the former is located 
closer to us, with the latter located furthe r away. Dis¬ 


tances of Cepheids show a bimodality (Mathewson et al. 
1986 Nidever et al.|2() f3j) th at is also observed in the HI 


velocities (Mathewson T984 Stanimirovic et al. 2004). 
Hence, each of our trends may be related to a different 
component. The trend at 65° is most likely present in 
the two components of the SMC. 

We obtained a turbulent component for the magnetic 
field of SB = (1.465 ± 0.069) /rG and for the or¬ 

dered magnetic field projected on the plane of the sky 
of B s k v = (0.947 ± 0.079) pG. The ordered-to-random 
field ratio at the SMC is closer to what is observed in our 
Galaxy than the average values for other irregular dwarf 
galaxies. 

This study is relevant for a better understanding of 
the magnetic field at the SMC, with a catalog contain¬ 
ing good polarization determinations, which can be used 
for several kinds of studies. We wish to emphasize that 
our data were concentrated at the NE and Wing sections 
of the SMC and part of the Magellanic Bridge. Further 
observations including the whole SMC, LMC, and Magel¬ 
lanic Bridge are necessary for a more complete picture of 
this system. The f orthcoming SOUTH POL project data 


(Magalhaes et ah 2012) will map the polarization for the 
sky. 


entire southern sky. It will cover the entire SMC-LMC 
system and Magellanic Stream and Bridge, which will al¬ 
low a better understanding of the spatial magnetic field 
behavior. SOUTH POL is expected to measure the po¬ 
larizations of around 24,500 bright stars (my < 15 mag) 
in this system with accuracy up to 0.1% and around 
218,000 stars (my < 17 mag) with accuracy up to 0.3%. 
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APPENDIX 

A. THE POLARIMETRIC CATALOGS 

Here we present a short version of the polarimetric cat¬ 
alogs. 

B. CLASSIFICATION OF THE FIELDS 

We used the F-test to quantify whether the number of 
Gaussians used to fit each of our fields was appropriate to 
represent the data set. The F-test applied for regression 
problems was used. We assessed whether a model with 
more parameters (one more Gaussian in our case) would 
fit the data better. The F statistic is given by 


( rss 1 -rss 2 



where RSSi is the residual sum of squares of model i. 
Pi is the number of parameters of model i, and n is the 
number of points used to fit the data. 

F has an F distribution with ( p 2 — pi, n—p 2 ) degrees 
of freedom. The null hypothesis of our test is that model 
2 does not fit the data better than model 1. We assumed 
a significance level probability of 2.5%, which gives sta¬ 
tistically significant results, and looked for the critical 
values of F in an online F distribution tablc]3 The null 
hypothesis is rejected if the F obtained from the data is 
greater than the critical value. Table [9] summarizes our 
results. 

The F-test demonstrates that all fields were well clas¬ 
sified. The negative values of F indicate that the RSS of 
model 2 is greater than that of model 1, which by itself 
already shows that model 1 fits the data better. The 
fields SMC 08, 10, 14, 25, 26, and 27 have limiting val¬ 
ues of F. In the case where RSS 2 » RSS\ such that 
RSS\/RSS 2 « 1, F simplifies to 

p = ip-p 2 ) / RSS 1 _ \ 

(P2 ~Pi)\ rss 2 > 

p _ (n-p 2 ) 

~ (P 2 -P 1 )' 

The six aforementioned fields satisfied this limit, indi¬ 
cating that the simpler model is the most robust. 

9 http://www.socr.ucla.edu/applets.dir/f_table.html 









































20 


Lobo Gomes, Magalhaes, Pereyra & Rodrigues 


Table 7 

Polarimetric catalog (observed polarization) 


ID 

R.A. 

(h:m:s) 

Decl. 

Pobs 

(%) 

a Pobs 

(%) 

PA 

(deg) 

V 

(mag) 

<T V 

(mag) 

0001 

1:00:05.29 

-71:49:46.23 

0.4120 

0.1330 

163.77 

18.58280 

0.16040 

0002 

1:00:04.48 

-71:52:41.15 

0.9270 

0.0490 

110.77 

16.71380 

0.16000 

0003 

1:00:03.63 

-71:54:54.25 

0.6460 

0.1050 

057.57 

17.59820 

0.16000 

0004 

1:00:05.28 

-71:51:45.71 

0.8860 

0.1710 

050.37 

17.31070 

0.16000 

0005 

1:00:05.76 

-71:50:21.05 

0.8300 

0.0930 

158.07 

18.65590 

0.16020 


Table 7 is published in its entirety in the electronic edition of the Astrophysical 
Journal. A portion is shown here for guidance regarding its form and content. 


Table 8 

Polarimetric catalog (intrinsic polarization) 


ID 

R.A. 

(h:m:s) 

Decl. 

(“:':") 

Pint 

(%) 

VPint 

(%) 

PA 

(deg) 

V 

(mag) 

(T V 

(mag) 

0001 

1:00:05.29 

-71:49:46.23 

0.5781 

0.1339 

179.78 

18.58280 

0.16040 

0002 

1:00:04.48 

-71:52:41.15 

0.6109 

0.0515 

110.33 

16.71380 

0.16000 

0003 

1:00:03.63 

-71:54:54.25 

0.8027 

0.1062 

046.57 

17.59820 

0.16000 

0004 

1:00:05.28 

-71:51:45.71 

1.0892 

0.1717 

043.28 

17.31070 

0.16000 

0005 

1:00:05.76 

-71:50:21.05 

0.9030 

0.0943 

168.31 

18.65590 

0.16020 


Table 8 is published in its entirety in the electronic edition of the Astrophysical 
Journal. A portion is shown here for guidance regarding its form and content. 


C. DEFINITION OF THE LIMITING 
MAGNITUDE FOR THE RANDOM 
BACKGROUND 

In order to define the limiting magnitude to filter the 
random background for the PFM fields, we made plots 
of the dispersion of 6 in function of the maximum my 
considered. We varied the limiting magnitude with an 
increment of 0.1 mag to construct the plots. The limit¬ 
ing magnitude was defined as the value where the curve 
linearly starts to increase, after having a saturation value 
(SMC12 and SMC20), or abruptly starts to increase 
(SMC13 and SMC23). Figure [II] shows these plots for 
the PFM fields. 
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Table 9 

F -test 


SMC 

Model 1 

Model 2 

77’ a 

-^data 

dof b 


Appropriate? 

01 

2 Gaussians 

3 Gaussians 

0.49 

(3,9) 

5.08 

yes 

02 

2 Gaussians 

3 Gaussians 

-1.17 

(3,9) 

5.08 

yes 

03 

1 Gaussian 

2 Gaussians 

-3.93 

(3,12) 

4.47 

yes 

04 

2 Gaussians 

3 Gaussians 

-1.34 

(3,9) 

5.08 

yes 

05 

Uniform 

1 Gaussian 

-1.85 

(2.14) 

4.86 

yes 

06 

1 Gaussian 

2 Gaussians 

-3.96 

(3,12) 

4.47 

yes 

07 

2 Gaussians 

3 Gaussians 

-2.71 

(3,9) 

5.08 

yes 

08 

2 Gaussians 

3 Gaussians 

-3.00 

(3,9) 

5.08 

yes 

09 

1 Gaussian 

2 Gaussians 

-1.77 

(3,12) 

4.47 

yes 

10 

2 Gaussians 

3 Gaussians 

-3.00 

(3,9) 

5.08 

yes 

11 

1 Gaussian 

2 Gaussians 

-1.21 

(3,12) 

4.47 

yes 

12 

1 Gaussian 

2 Gaussians 

-2.74 

(3,9) 

5.08 

yes 

13 

1 Gaussian 

2 Gaussians 

-2.16 

(3,11) 

4.63 

yes 

14 

Uniform 

1 Gaussian 

-7.50 

(2,15) 

4.76 

yes 

15 

1 Gaussian 

2 Gaussians 

-1.88 

(3,12) 

4.47 

yes 

16 

2 Gaussians 

3 Gaussians 

-0.87 

(3,6) 

6.60 

yes 

17 

2 Gaussians 

3 Gaussians 

-1.88 

(3,6) 

6.60 

yes 

18 

1 Gaussian 

2 Gaussians 

3.00 

(3,12) 

4.47 

yes 

19 

Uniform 

1 Gaussian 

-0.07 

(2.14) 

4.86 

yes 

20 

1 Gaussian 

2 Gaussians 

-1.58 

(3,10) 

4.83 

yes 

21 

1 Gaussian 

2 Gaussians 

-0.22 

(3,9) 

5.08 

yes 

22 

2 Gaussians 

3 Gaussians 

-1.99 

(3,8) 

5.42 

yes 

23 

1 Gaussian 

2 Gaussians 

-1.98 

(3,7) 

5.89 

yes 

24 

1 Gaussian 

2 Gaussians 

-2.65 

(3,8) 

5.42 

yes 

25 

Uniform 

1 Gaussian 

-6.50 

(2,13) 

4.96 

yes 

26 

1 Gaussian 

2 Gaussians 

-2.33 

(3,7) 

5.89 

yes 

27 

2 Gaussians 

3 Gaussians 

-2.33 

(3,7) 

5.89 

yes 

28 

1 Gaussian 

2 Gaussians 

-1.26 

(3,12) 

4.47 

yes 


a F value obtained from the data. 
b Degrees of freedom (p 2 — pi, n — P 2 ). 
c F critical value. 
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